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ABSTRACT 

Large amplitude, electrostatic plasma waves are relevant to physical 
processes occurring in astrophysical magnetospheres wherein charged 
particles are accelerated to relativistic energies by strong waves 
emitted by pulsars, quasars or radio-galaxies. In this paper the non- 
linear, relativistic theory of traveling Langmuir waves in a cold 
plasma is reviewed. The cases of streaming electron plasma, electron- 
ion plasma and two.-streams are discussed. 


1. Introduction 

Strong waves capable of driving plasma particles to relativistic 
energies are of current interest in several fields of research. A 
wave is strong enough to drive electrons or positrons relativistic if 
its dimensionless strength parameter v = eE ^ Y /m p ac is comparable to 
unity. " 

In astrophysical context, it was first proposed by Pacini*^ that 
a rotating neutron star with misaligned magnetic dipole and rotation- 
al axes can release its magnetic and rotational energy by emitting 

low-frequency magnetic dipole electromagnetic waves at the stellar 
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rotation frequency. Gunn and Ostriker J showed that charged par- 
ticles could be accelerated to ultrarelativistic energies by low- 
frequency pulsar magnetic dipole radiations with v > 10 1 , thus sug- 
gesting that pulsar electromagnetic fields may be a cosmic ray source. 
Kennel et al 1 found that when ions as' well as electrons are driven 
relativistic by a plane plasma wave, there is an upper limit to the 

cosmic ray number flux above which the wave encounters cut-off. 
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Rees suggested that the rotational energy of a large number of pul- 
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sars or collapsing stars could supply the power for quasars or radio- 
galaxies through the emission of low-frequency strong electromagnetic 
waves. 

Recently, Bodo et al S ^ used the dispersion relation of strong elec 
tromagnetic waves to study the propagation conditions of high-frequen 
cy coherent radio emission in pulsar magnetospheres. Chian and 
Kennel^ estimated that the strength parameter of coherent pulsar ra- 
diowaves in the source region lies in the range 10 -2 <v<l0 3 ; in the 
field of those waves the charged particles can acquire from weakly to 
moderately relativistic energies. 

It is worth mentioning that in addition to astrophysical applica- 
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tions, strong waves can be relevant to laser-plasma interaction ’ 
and particle acceleration by laser beatwave * . At laser intensities 
of 10 16 Wcm -2 for CO 2 gas laser and 10 18 Wcm -2 for Rd glass laser the 
wave strength parameter v - 1 and relativistic effects become impor- 
tant. Lasers of such intensities are already available in laborator- 
ies. 

Astrophysical applications of strong electromagnetic waves have 

been reviewed elsewhere"^ * . In this paper only the subject of 

nonlinear, relativistic electrostatic Langmuir waves is reviewed. 

The Langmuir mode may play an important role in both emission and 

propagation processes of pulsar radiations. 
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Ruderman and Sutherland suggested that the two-stream instabil- 
ity resulting from the relative streaming between electrons and posi- 
trons in a pulsar magnetosphere can excite unstable electrostatic 
plasma waves which then bunch the charged particles, producing coher- 
ent microwave radiation. Khakimova and Tsytovich^^ and Buti^^ con- 
sidered the generation of coherent pulsar radio emission by Langmuir 
solitons arising from the modulational instability that develops 

from electrostatic plasma waves if their energy exceeds a certain 
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threshold. Mamradze et al 1 showed, that the coalescence of two Lang 
muir waves can lead to infrared pulsar emission. Mikhailovskii^^ 
presented a theory of pulsar radio emission using the decay process 
of a Langmuir wave into two electromagnetic waves. 

As a strong electromagnetic wave propagates in a pulsar magneto- 
sphere, its absorption and scattering is often dominated by collec- 
tive processes involving electrostatic plasma waves. Tsintsadze and 



Tsikarishvili 17 ^ investigated the parametric instability of a strong 
electromagnetic wave decaying into Langmuir and ion-acoustic waves. 
Drake etal. 18 ^ suggested that such instability may break up strong 
electromagnetic waves before energy transfer to magnetospheric par-' 
ticles 2 ^’ 3 * * * * * ^ can occur. Arons et al.* 9 '* considered the parametric 
Raman scattering of strong pulsar radiations by subluminous Langmuir 
waves. 

The subject of Langmuir waves has been studied extensively. To the 
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linear approximation, Tonks and Langmuir ' demonstrated that a small 
displacement of a uniform plasma results in an electrostatic force 
which gives rise to harmonic oscillations at the plasma frequency. 

As the amplitude of oscillations increases and the relativistic vari- 
ation of particle mass appears, many nonlinear effects can alter sig- 
nificantly the behavior of Langmuir waves. As an initial step to un- 
derstand the complex role that the Langmuir mode plays in various 
strong wave problems, we examine the traveling wave solutions of non- 
linear, relativistic Langmuir waves. Three cases are studied: (i) 

streaming electron plasma, (ii) electron-ion plasma and (iii) two- 
streams. 

2. Streaming Electron Plasma 

The traveling wave theory of nonlinear, relativistic, Langmuir 
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waves was first studied by Akhiezer and Polovin J . The case they 

treated was a stationary electron plasma. In many astrophysical prob- 

lems one encounters electron (or positron) streams moving at relativ- 

istic stream velocities. In this section the treatment of Akhiezer 
and Polovin 2 "^ is generalized to the case of streaming electron 
plasma . 

Consider a large-amplitude Langmuir wave propagating with a con- 
stant speed V in a cold, unmagnetized plasma. The average electron 
charge and current densities are neutralized by positive ions, but 
the ion dynamics is ignored. Since we are seeking traveling wave 
solutions we can write the ba.sic equations in terms of a combined 
space-time variable 0 = t - nx/c, where the index of refraction 
n(-<» < n < «) is a- measure of the wave speed V = c/n. The basic equa- 
tions are the relativistic equation of motion, the continuity equa- 
tion and Poisson's equation 
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supplemented by Maxwell's equation 
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where y = (1 - v 2 /c 2 ) ^ and J q are neutralizing ion charge and 
current densities. In nonlinear theory, the plasma stream velocity 
can be defined as 

V = <Nv>/ <N> , ( 5 ) 


which can be obtained by taking the phase-average of (3) and (4), 

with J = eN V . 
o os 

The basic equations (1) - (5) can be combined to yield a single 
equation 


i.du, 2 n - y + y/c 

2 dt^ (1 - nu/y) 2 


where u = yv/c = reduced electron velocity, x = w 0, <d 2 = 4irN e 2 /m 

pe pe o 

and W is a constant that characterizes the wave amplitude. Evidently 
(6) indicates that periodic wave solutions exist. The nonlinear dis- 
persion relation is obtained by a further integration of (6) 
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where P is the wave period. (7) describes Langmuir wave oscillations 
for arbitrary wave amplitudes. It recovers the result of a station- 
ary plasma if V s is set to zero, and can be obtained from the dis- 
persion relation for a stationary plasma by a Lorentz transforma- 
22 ) 

tion . Fig. 1 shows the variation of the cut-off frequency (n = 0) 


of nonlinear, relativistic Langmuir waves with the amplitude of elec- 
tron oscillations (where v^ = |vj - V 2 I, v 1>2 being the turning points 
of v).for different stream velocities. It demonstrates that, due to 
relativistic variation of particle mass (m ■+■ ym ) , the Langmuir wave 
frequency decreases as the wave amplitude increases. In Fig. 2, the 

product of c 2 /u 2 with k 2 (where k = nrn/c) is plotted against the 
P e , 

wave speed for different wave amplitudes. It indicates that, for a 
given wave velocity, relativistic effects cause the wavenumber of 
Langmuir waves to decrease as the wave amplitude increases. It also 
shows that there are regions in the dispersion curves where periodic 
wave oscillations are not possible due to the phenomenon of wave- 
breaking. 

The wave speed of nonlinear, relativistic Langmuir waves can be 
either superluminous (V > c) or sub luminous (V < c) . guper luminous 
waves can attain arbitrary amplitudes, but the amplitude of sublumi- 
nous waves is limited. . For nonlinear subluminous waves moving at a 
given speed, there is a critical wave amplitude above which wave 
breaks whereby ordered wave energy is transformed to random plasma 
thermal energy. In laser fusion research it has been found that 
wavebreaking of electrostatic plasma waves can produce highly ener- 
getic electrons 2 **^ ; similar processes are expected to take place in 
astrophysical plasmas. The problem Of wavebreaking can be clarified 
by analyzing the first integral corresponding to (6) in the time-inde- 
pendent frame 22 ^’ 2 ^ 


1 dV 2 

10 ="-f(Y) 


( 8 ) 


with 


fi Y) = | 1/n - V s /c|(y 2 - l)* 5 + (1 - V s /nc)Y , (9) 

and t = (1 - n . z/c. The mathematical form of (8) suggests 

that the nonlinear wave problem is analogous to the motion of a par- 
ticle of unit mass at position y, at time t, in Newtonian dynamics. 
Thus, (8) corresponds to the conservation of energy with (dY/dr) 2 /2 = 
kinetic energy, f(y") = potential energy and W = total energy. A 

4 

typical plot of / as a function of y is displayed in Fig. 3, which 


shows that there is a maximum value of W for a periodic solution to 
exist 


V = f(X = 1) = 1 - V / nc , (10) 

max s 

above which wavebreaking then occurs. Transforming the wavebreaking 
condition back t'o the laboratory frame, it shows that traveling Lang- 
muir wave solutions exist only if the electron velocity lies in the 
range 



v < v < c/n 

if V < V 



< 

c .. 

c/n < v < v 
L- c 

s 

if V > V 
s 

i 

(ii) 


where 

1 + V’ 2 /c 2 - 2nV /c 

v /c = ^ 5 . (12) 

c 2V /c - n(l V 2 /c 2 ) 

s s 

Alternatively, the wavebreaking condition can be expressed in terms of 

R 2 = (eE /mew ) 2 as 
v max pe 

R 2 > R c 2 = 2[(1 - n~ Z y h (1 - V s /nc) - (1 - V^/c 2 )* 5 ] . (13) 

Graphical displays of (11) and (13) axe given by Fig. 4, where the 
shaded areas represent the regions where wavebreaking occurs. We see 
that subluminous Langmuir waves can attain large amplitudes if their 
speeds are close to the speed of light; as the wave speed gets away 
from c the critical wave amplitude decreases and becomes zero if 



Typical waveforms are shown in Fig. 5 for the = 0 case. It is 
seen that, as v^ approaches the critical value c/n the waveforms 
steepen, indicating the possibility of wavebreaking. For low wave 
amplitudes (v < 1) waveforms are essentially sinusoidal, whereas for 
large wave amplitudes (v » 1) waveforms become triangular. 

3. Electron- Ion Plasma 

In the previous section, only the electron motion was considered 


whereas the motion of the positive ions was ignored. The results ob- 
tained thus apply to comparatively small wave amplitudes for which the 
ion dynamics are insignificant because of its large rest mass. In the 
presence of very large-amplitude waves, ions can attain considerable 
directed velocities, becoming relativistic if yv » 1, where y = m /nn 
Extremely intense sources of electromagnetic radiation do exist in 
astrophysical plasmas. For instance, in a Crab pulsar magnetosphere, 
the strength parameter v associated with the low-frequency magnetic 

dipole radiation can be > 10 11 and the ions can acquire ultrarelativ- 
2 ) 

istic velocities 

Consider now a two-fluid electron-ion plasma in the absence of ex- 
ternal fields 2 ^ ~ 2 ^ . In this case the basic equations for Langmuir 
waves are the relativistic equations of motion and of continuity for 
electrons and ions, respectively, and Maxwell's equations 


- = 4ire(N. - N ) 

c d0 v 1 e' 


+ 4re(N.v. - N v ) 
d6 i i e e J 


= 0 


(14) 

(15) 


Taking the phase-average of (14) and (15), we see that in the absence 
of external fields the stream velocities of electrons and ions are 
equal 


V = <N v >/<N > = <N.v.>/<N.> . (16) 

s ee e l l l 


Stationary electron-ion plasma corresponds to the particular case in 
which the average particle fluxes are zero. 

The basic equations can be combined to give two conservation rela- 
tions, one expresses conservation of energy 2 ^ 


E 2 /8ir + Nm y c 2 + Nm.y.c 2 = Nm c 2 W , 

e e . li e 

and another expresses conservation of particle momentum 


(17) 


(u. + yu ) + y (u_ - ny ) = D , (18) 

A C C C 

where N = <N>(1 - nV s /c) , u =- yv/c, W and D are constants. The dis- 


persion relation for nonlinear, relativistic, Langmuir waves in an 
electron-ion plasma can be obtained from (17) and (18) 


I 



,u 2 1 - nu e /T e 
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(W - Y 

1 uj v e 
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When evaluating a particular solution for an electron-ion plasma 
it is necessary to decide what values should be assigned to the two 
conservation constants W and D. The constant W, as shown in the pre- 
vious section, can be considered as a wave amplitude parameter. How- 
ever, the role of the constant D is more subtle. Most papers treat D 
either as an arbitrary constant or as a plasma drift velocity. This 
leads to mathematically correct, nonetheless unphysical solutions. In 
the absence of a wave, v = v. = V , (18) gives 
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(l+-y)(l - V s 2 /c 2 ) ** V.C- n(l + y) (1 - V^/c 2 )" 1 * = D , (20) 


which shows that'D is a function of V^, n and y. For large wave amp- 
litudes', D also becomes a function of the wave intensity. In fact, it 
was established that D is uniquely related to the amplitude parameter 
w 30 ). Hence, in order to obtain a specific solution with given plasma 
stream velocity, wave velocity, electron-ion mass ratio and wave in- 
tensity, the dynamic parameter D must be adjusted accordingly so as 
to render self-consistency to the solution. In Fig. 6, the variation 
of the self-consistent dynamic parameter D with u (where u = |ui~U 2 |) 
for a stationary electron-ion plasma is displayed for different 
values of n. 

The general behavior of nonlinear, relativistic Langmuir waves in 
an electron-ion plasma is similar to the case of an electron plasma. 
The major effect of ion dynamics is to decrease the period of oscil- 
lation; this effect gets more pronounced as the wave amplitude in- 
creases. In fact. Kennel and Pellat^^ showed that for extremely 
large wave amplitudes the contribution of ions toward wave properties 
is equal to electrons; when the wave energy density greatly exceeds 
the rest mass energy, the inertia of plasma particles is not deter- 
mined by their rest mass but by the kinetic energy they acquire from 
the waves. Hence, in the limit of extremely large wave intensities, 
an electron-ion plasma behaves like an electron-positron plasma. 


4. Two-Streams 

Two-stream instability can be of fundamental importance to the gen- 

12 ) 

eration of coherent pulsar - radio emission . In this section', we 
briefly discuss the relativistic effects in the two-stream instability 
by analyzing the traveling wave solution of nonlinear, relativistic 
Langmuir waves. As an illustration, we treat the case of two inter- 
penetrating electron and ion streams. In order to permit relative 
streaming between electrons and -ions, (14) and (15) are modified to- 


ll ~ 
c d6 


4tt (eN ^ 


eN + p ) 
' e o J 
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de 
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eN v + J ) = 0 
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where p and J are respectively the constant charge and current den-- 
0 o 

sities required to neutralize the average charge and current densities 
of the streams. Evidently, (21) and (22) show that electron stream 
velocity V is different from ion stream velocity V . , namely 

SC SI 


<N v >/<N > ^ <N. v. >/<N. > 
-- e e e - i i i 


with 


(23) 


p = e(<N > - <N. >) , J = e (<N v > - <N.v.>) . (24) 

o v e l o v ee 11 

The dispersion relation for this case can be obtained by making 

22) 25) 

use of the Lorentz transformation technique . The detailed re- 

31) 

suits will be presented in a forthcoming paper . Fig. 7 shows some 
examples of dispersion curves for electron-ion streams, with various 
wave amplitudes. We see that the two-stream instability occurs for •. 
wavenumbers below a certain critical value. Similar to the case of 
a single stream (see Fig. 2), for a given wave velocity, relativistic 
effects cause the wavenumber to decrease as the wave amplitude in- 
creases. Hence, for nonlinear, relativistic Langmuir waves in 
electron-ion streams, there is a range of wavenumbers that are un- 
stable according to the linear theory- but becomes stable as the wave 
amplitude gets sufficiently large. Moreover, Fig. 7 indicates that 
wavebreaking plays an important role in determining the dispersion 


characteristics of electron-ion streams. For, as the wave amplitude 
increases, the range of wave velocities for which wavebreaking occurs 
also increases. If the wave amplitude becomes very large, wavebreak- 
ing may occur for all subluminous wave velocities that lie between the 
electron and ion stream velocities. 


5. Conclusion 

The traveling wave solutions of nonlinear, relativistic, Langmuir 
waves studied in this paper serve to understand some basic properties 
of large-amplitude electrostatic plasma waves when v/ave intensity in- 
duced relativistic and nonlinear effects are important. Astrophysical 
implications of these wave properties for strong wave problems remain 
to be explored. In addition, many effects on nonlinear, relativistic, 
Langmuir waves such as .plasma temperature*^ and wave stability 

35), 36)^ left out in this paper, await further investigations. 
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Figure 1. Variation of u/u) with v /c for indicated values of V /c 
n P e a s 
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V/c 


Figure 2. Variation of (c 2 /w 2 )k 2 with V/c for indicated values of 

v/c; V = 0.5c. pe 
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Figure 5. Variation of E/E max with t for (a) v Q = 0.5c(u o = 0.58) 
(b) v 0 = 0 . 995c(uq =10). 
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